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The realization of p—n junctions in graphene, combined with the gapless and chiral nature of its 
massless Dirac fermions has led to the observation of many intriguing phenomena such as quantum 
Hall effect in bipolar regime, Klein tunneling, and Fabry-Perot interferences all of which involve 
electronic transport across p—n junctions. Ballistic snake states propagating along the p—n junctions 
have been predicted to induce conductance oscillations, manifesting their twisting nature. However, 
transport studies along p—n junctions have so far only been performed in low mobility devices. 
Here, we report the observation of conductance oscillations due to ballistic snake states along a p-n 
interface in high quality graphene encapsulated by hexagonal boron nitride. These snake states are 
exceptionally robust as they can propagate over 12 pm, limited only by the size of our sample, and 
survive up to at least 120 K. The ability to guide carriers over a long distance provide a crucial 
building block for graphene-based electron optics. 


Snake states, a counterpart of skipping orbit states, 
were first proposed to exist along the boundary be¬ 
tween positive and negative magnetic field [TH5]. How¬ 
ever, the difficulty to achieve a magnetic field gradient 
experimentally [6] has led to a new proposal based on a 
p-n junction in graphene J. 8] where, instead of inverting 
the magnetic field polarity, we interconvert carrier types 
between hole and electron across the junction, generat¬ 
ing snake-like trajectories at finite magnetic field. An in¬ 
crease in conductance along the p-n interface of graphene 
devices on Si 02 at finite magnetic field has been observed 
and attributed to snake states [SJ. It has been theoret¬ 
ically predicted that, in a higher quality device, these 
snake states, even though propagating along the p-n in¬ 
terface, can manifest themselves through oscillatory con¬ 
ductance across the interface as a result of alternating 
end points of the snake states [71 H0HI3] . The high-quality 
p-n junction in graphene, only feasible recently follow¬ 
ing the development of graphene heterostructures fT4l H5] . 
provides a perfect platform for snake states due to its 
gapless and chiral nature which leads to a transparent 
p-n junction at normal incident angle, owing to Klein 
tunneling [lBHIS] . 

In this work, we investigate conductance oscillations 
arising from ballistic snake states by employing a top¬ 
gated graphene encapsulated by hexagonal boron nitride 
(hBN) in order to access the ballistic regime |T4, 15, [20!. 
The observed conductance oscillations are present only at 
finite magnetic field and in bipolar regime. The density 
and magnetic field dependence of these oscillations sup¬ 
ports a theory of snake states propagating along the p-n 
junction. The oscillations persist up to at least 120 K, 
showing the robustness of the snake states. 

Figures la and b show a schematic diagram and an 


atomic force microscopy (AFM) image of our device. 
Graphene is transferred onto a 32-nm-thick hBN and 
etched into a Hall bar geometry whose width and length 
is 900 nm and 12.2 /im respectively. The device is then 
encapsulated by an 8-nm-thick hBN. A top gate is de¬ 
fined onto the device, parallel to the Hall bar (see Meth¬ 
ods for fabrication). The thin top hBN flake allows us 
to create a sharp p-n interface which increases the trans¬ 
mission probability across the p-n junction [17 . We mea¬ 
sure conductance (G) using a two probe method in which 
carriers are injected into graphene through an injector (i) 
and collected at a collector (c). Back gate (BG) and top 
gate (TG) allow us to control densities in back-gated and 
top-gated regions (BGR and TGR) independently. 

Under a finite magnetic field normal to graphene sur¬ 
face (B = Bz), carriers, injected into graphene, will un¬ 
dergo a skipping orbit motion along the edge (Fig. lc, 
blue curve). After reaching a p-n interface, some of 
the carriers will convert into snake states (Fig. lc, black 
curve), propagating along the interface and inducing con¬ 
ductance oscillations. Figure Id shows conductance as a 
function of top gate voltage (Utg) at Vbg = 17.2 V, when 
the BGR is n-doped. As we sweep Vtg, carrier density 
in the TGR can be tuned from p-type doping to n-type 
doping, creating n-p or n-n' junctions respectively. At 
B — 0 T (black line), the conductance shows no oscilla¬ 
tions as a function of Vtg in both unipolar and bipolar 
regime. However, at B = 2 T (red line), clear oscillations 
emerge in the bipolar regime (—4 V < Vtg < -1 V, 
thick red line) but is absent in the unipolar regime 
(1 V < Vtg < 4 V). When BGR is p-doped (Fig. le), 
similar oscillations can also be observed at finite field in 
the bipolar regime (1 V < Vtg < 4 V). 

The oscillations we observe here only exist in the bipo- 
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FIG. 1: Device schematic and conductance measurements, (a) Schematic diagram of a graphene heterojunction device. 
Graphene is encapsulated by hBN and top gate and back gate are used to create a p-n junction, (b) False-colour AFM image 
of a device after mechanical cleaning with an AFM tip. The contacts labeled by i and c are used for an injector and collector 
with the width of ~300 nm. The width and length of graphene are 900 nm and 12.2 pm respectively. A top gate (orange) 
is overlaid whose width is 500 nm. The edge of the top gate is aligned with the graphene edge. The scale bar is 2 pm. (c) 
Classical trajectories of skipping orbits (blue) and snake states (black). Magnetic field B = B z is applied normal to graphene 
surface where z is the unit vector normal to the surface, (d, e) Conductance versus Vtg at Vbg = 17.2, —20 V respectively. 
Black and red curves represent conductances at 0 and 2 T respectively. Conductance oscillations (thick red line) can be seen 
only at finite field in the bipolar regime. Curves are offset for clarity. 


lar regime which allows us to rule out Shubnikov-de 
Haas (SdH) oscillations as a main mechanism. Similar 
oscillations in the bipolar regime are also observed in 
graphene heteroj unctions as a result of Fabry-Perot (FP) 
interferences[19 . However, one distinct difference is that 
our oscillations are absent at zero magnetic field while 
FP oscillations can be observed at zero field. This dis¬ 
misses the possibility that our oscillations originate from 
the FP interferences. 

In order to gain more insight into the nature of these 
oscillations, we measure conductance as a function of Vbg 
and Vtg- To eliminate the background conductance, the 
numerical derivative of G with respect to Vtg is ana¬ 
lyzed (see Supplementary Fig. 1 for data before numer¬ 
ical derivatives). At B = 0 T (Fig. 2a and b), we ob¬ 
serve two lines which divide (Vbg, Vtg) plane into four 
quadrants. The vertical line corresponds to the density 
in the BGR (tibg) equal to zero while the diagonal line 
follows the zero density in the TGR (titg)- The cross¬ 
ing between these two lines at (V^g, ^tg) = ( — 3.2, 0) V 
indicates the charge neutrality point of the full sample. 
The four quadrants correspond to rm', pn, pp', and np 
doping where the first and second letters are carrier types 
of the BGR and TGR respectively. Evidently, no oscil¬ 
lation is present at zero field. As we increase B to 2 T 
(Fig. 2c and d), oscillations become apparent in the pn 
and np quadrants. The curvatures of these oscillations 


are distinctively different from FP oscillations jT9j or SdH 
oscillations in the presence of p-n iunctions [2TH23] . 

The combination of magnetic field and a p-n junc¬ 
tion leads us to consider the contribution from the snake 
states. Figure 2e shows a diagram of a snake state prop¬ 
agating along the p-n interface. Here, we assume a step¬ 
like potential at the interface. With applied magnetic 
field, the carrier trajectory is bent into a circular orbit 
with cyclotron radius r c = Hk^/eB = hy/wn/eB where 
h is the reduced Planck’s constant, kp is the Fermi wave 
vector, n is the carrier density, and e is the elementary 
charge. The oscillations in conductance can be under¬ 
stood by considering the end point at which a snake state 
exits the p-n junction. For instance, if a snake state ends 
up on the collector side (Fig. 2e), this will increase the 
transmission probability from the injector to the collec¬ 
tor, enhancing the conductance. However, if it exits the 
p-n interface on the injector side (Fig. 2f), the conduc¬ 
tance should diminish. The choice of the end point of the 
snake states will depend on the starting point, the length 
of the p-n interface (L) and the cyclotron radii of BGR 
(rf G ) and TGR (rJ G ) which vary with B and n. For 
example, we can move the end point of the snake state 
in Fig. 2e to the one in Fig. 2f by increasing titg- This 
enlarges the cyclotron radius in the TGR and causes the 
snake state to end up on the injector side instead of the 
collector side. 
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In the bipolar regime, we can determine the end point 
of a snake state by considering the following conditions. 
For Lo = 2N(rf G + rJ G ) + rf G where N is the integer 
part of the ratio L/ (2r^ G + 2rJ G ), a snake state will end 
on the collector side if L > Lq and on the injector side 
if L < Lq. Using these conditions, we simulate the end¬ 
point in the bipolar regime as a function of Ubg and Utg 
where 77 -bg = CbgV-qq, tttg = GbgFbg + GtgUtq, G 
is the capacitive coupling, and U* = V — U D . Figure 2g 
shows the result from such simulation at B = 2 T. Here, 
we make the assumption, based on the semiclassical bil¬ 
liard model with a finite size injector[12], that all snake 
states have the same starting point. We note that there 
is no fitting parameter in our simulation as all parame¬ 
ters (L, I>, tibg and utg) can be determined experimen¬ 
tally. The simulation agrees relatively well with the data, 
especially at high density. The discrepancy at low den¬ 
sity is likely due to disorder as cyclotron radius becomes 
smaller and comparable to the width of the depletion 
region of the p-n junction (see below for detailed discus¬ 
sion). The contour line of a constant dG/dUxG is there¬ 
fore described by r;? G +rJ G oc yj tibgH - \J ^tg = constant, 
explaining the curvature of the oscillations observed in 
Fig. 2g. We note that, under the experiment conditions 
we apply, the snake states need to transmit at least ^50 
times across the p-n interface in order to propagate from 
one end to the other end which is about 12.2 pm long. 
This shows the robustness of the snake states compared 
to the skipping orbit states along a graphene edge in 
which the focusing peaks of the transverse magnetic fo¬ 
cusing are greatly diminished after a few reflections of the 
edge [2m [24] . While the skipping orbits along the edge of 
graphene have to suffer potential fluctuation from dan¬ 
gling bonds, a lack of dangling bonds as well as a colli¬ 
mating effect of the p-n junctions possibly contribute to 
the stability of the snake states. 

In the unipolar regime, the absence of the oscillations 
can also be explained by similar consideration of the car¬ 
rier trajectories along the interface. Since both sides of 
the interface have the same type of carriers, they circulate 
in the same direction. However, the relative magnitude 
of the densities on both sides produces two distinct types 
of trajectories. For tibg > ^tg (Fig- 2h), we have the cy¬ 
clotron radius in the BGR larger than that in the TGR 
and the carrier can drift along the p-n interface simi¬ 
lar to the snake states. However, the main difference is 
that carriers in this density configuration always end up 
on the injector side. Hence, the transmission probability 
remains unchanged as Ubg and Vtg are varied. When 
nBG < ^tg (Fig. 2i), the carrier will transmit through 
the p-n interface without drifting along the interface and 
always ends up on the collector side. As a result, no os¬ 
cillation is achieved in both cases. 

We now turn to the magnetic field dependence of the 
oscillating conductance. Figure 3a shows dG/dUxG as a 
function of Utg and B at Ubg = —20 V. The conductance 


oscillations in the bipolar regime (0.5 V < Utg < 4 V) 
are apparent from T onward and they are much more 
well-behaved at low field (< 2 T) than at high field. To 
understand the observed feature, we simulate the end 
points of the snake states based on the aforementioned 
conditions. The simulation illustrated in Fig. 3b shows 
that the oscillations are periodic in B for a fixed Utg 
with the period given by 12 

2h 

A B = — (v^BG + v^tg) (1) 

The simulation matches the data reasonably well at 
low field. However, it fails to capture the onset of the 
conductance oscillations and aperiodicity at high field. 
The onset seems to relate to the snake states undergo¬ 
ing a deconfinement transition below a critical field [25] 
B c = hy/wn/eW where W is half the width of the TGR. 
As the magnetic field is lowered, the cyclotron radius 
becomes larger. Consequently, once the radius is com¬ 
parable to half the width of the TGR, the snake states 
can back scatter to the other side of the TGR edge caus¬ 
ing the oscillations to vanish. The fit of B c to the onset 
of our data yields (blue dashed lines in Fig. 3a and b) 
W = 180 nm which is roughly half the width of the TGR 
(^250 nm). Note that a similar argument can also be 
applied to the BGR which yields B c = 0.57 T for the 
200-nm-wide half width (vertical line in Fig 3b). More¬ 
over, the onset as a function of density at B = 1.5 T 
agrees quite well with the line of constant rJ G = 180 nm 
(see Supplementary Fig. 2 and Supplementary Note 1). 

The aperiodicity at high field is likely related to dis¬ 
order. As the field is increased, carriers propagate closer 
to the depletion region of the p-n junction. In this re¬ 
gion, the effect of the disorders becomes prominent be¬ 
cause of the reduced screening at low carrier density [26.. 
The density fluctuation 127, [28] increases the roughness 
of the p-n interface which leads to complex carrier 
trajectories[29, [30] . The interface roughness can reflect 
the incoming snake states due to exponential suppression 
of transmission rate at finite incident angle [17]. There¬ 
fore, it acts as a scattering site which behaves as a new 
injector on the p-n interface, establishing a shorter effec¬ 
tive length (L e ff) the snake states needed to travel. The 
inset in Fig. 3b shows L e s calculated from the period in 
Fig. 3a (the difference in B between two adjacent min¬ 
ima) at Utg = 4 V. At low field, the cyclotron radius is 
still large and most snake states can skip over scattering 
sites, resulting in L e ^ comparable to the length of our 
p-n interface. As B is increased, the snake states now 
with smaller cyclotron radius cannot avoid the scatter¬ 
ing sites and, as a result, L e ^ decreases. The simulations 
in which we model the effect of disorder by putting a 
scattering center of finite size on the p-n interface show 
similar increasing periods with increasing B as observed 
in the data (see Supplementary Fig. 3 and Supplemen¬ 
tary Note 2). 
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FIG. 2: Density dependence of the oscillating conductance and carrier trajectories along the p—n interface, (a) 

dG/dVrG as a function of Vtg and Vbg at 0 T. No oscillation is observed, (b) Vertical cuts in a at Vbg = 17.2 V (red) and 
—20 V (black), (c) dG/dVrG as a function of Vtg and Vbg at 2 T. Conductance oscillations are clearly present in the bipolar 
regime, (d) Vertical cuts in c at Vbg = 17.2 V (red) and —20 V (black). The red curves in b and d are offset for clarity, 
(e, f) Classical trajectories in the pn density configuration. By tuning carrier densities, we can vary a cyclotron radius and 
alternate the end points of the snake states between the injector and collector sides, (g) Simulation of the end points of snake 
state trajectories along the p-n interface as a function of Vtg and Vbg- When a trajectory ends on the collector (injector) side, 
the trajectory is represented in orange (blue), (h, i) The carrier trajectories in the nn' density configuration. In h (i), the 
trajectories always end on an injector (collector) side. 


Finally, we comment on the temperature dependence 
of the conductance oscillations. Figures 3c-e display the 
conductance oscillations in the pn regime at 5, 60, and 
120 K respectively. The presence of the oscillations at 
such elevated temperature supports our claim that, in¬ 
stead of quantum interference, they are induced by snake 
states, which rely on ballistic transport and appearance 
of caustics[TT]. Simple estimation based on temperature 
broadening of k F m eb yields that quantum interfer¬ 
ence should be visible only up to ~10 K. Some of the 
irregular features at low temperature such as phase shift 
(Fig. 3c, red arrow), or bifurcations (Fig. 3c, red cir¬ 
cle) also disappear at high temperature. This suggests 
that such features could be related to quantum interfer¬ 
ence between different snake states [32]. More theoretical 
study is needed to understand the origin of these fea¬ 
tures. In addition, a new fabrication method based on 
the van der Waals assembly process [33], which produces 
a large area of encapsulated graphene without bubbles 


and wrinkles, could potentially lower the onset of the 
conductance oscillations in B and preserve the ballistic 
motion up to room temperature. This may lead to the 
room temperature application of the snake states such 
as electron waveguides [4] and electrically reconfigurable 
wiring [29]. 

The observation of the conductance oscillations is quite 
surprising given that the size of our injector is larger than 
r c which could have led carriers to exit the interface on 
both sides equally. However, under magnetic field, car¬ 
riers injected into graphene will undergo a skipping or¬ 
bit motion which restricts them to propagate along the 
edge [2(1. The lateral extension of the skipping orbit mo¬ 
tions will be on the order of r c which can be much smaller 
than the size of our injector (^300 nm). In addition, 
due to the exponential collimation of a graphene p-n 
junction 17, [19], only carriers which encounter the p-n 
interface near normal incident angle will be able to pass 
through, owing to Klein tunneling. For carriers which 
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FIG. 3: Magnetic field and temperature dependence of the oscillating conductance, (a) dG/dVrG as a function 
of Vtg and B at Vbg = — 20 V. The conductance oscillations are regular at low field and become more chaotic at high field. 
The blue dashed line is the critical magnetic field in the TGR. (b) Simulation of the end points of snake state trajectories as a 
function of Vtg and B at Vbg = —20 V. The shaded regions represent the area in which magnetic field is below the critical field 
in the BGR (vertical slab) or TGR. Inset: The effective lengths of the p-n interface versus B at Vbg = — 20 V and Vtg = 4 V. 
Error bars, calculated from Eq. 1, represent the uncertainty in determining the period of magnetic field, (c, d, e) dG/dVrG in 
the pn density configuration at 5, 60, and 120 K respectively. The arrow and circle in c indicate phase shift along the diagonal 
line and bifurcations respectively. The colour-scale bar is applied to all c, d, and e. 


collide the interface at large incident angle, they will 
mostly be reflected back and behave as a skipping or¬ 
bit states along the p-n interface. Therefore, the effect 
of a finite size injector is likely to introduce large conduc¬ 
tance background and reduce amplitude of the conduc¬ 
tance oscillations as only a small portion of electrons will 
turn into snake states. Moreover, simulations based on 
the semi-classical billiard model of a graphene p-n junc¬ 
tion have shown that the conductance oscillations still 
persist even for an injector whose width is wider than a 
cyclotron radius l [T2 ; T3J. 

Furthermore, we observe that the conductance oscilla¬ 
tion pattern is independent of an injector location (see 
Supplementary Fig. 4). A wide injector acts as a line 
source, instead of a point source, producing all possible 
skipping orbits. Hence, the distribution of electron tra¬ 
jectories before encountering a p-n junction should be 
identical along graphene edge. This leads to the same 
initial conditions for incident electrons at the p-n junc¬ 
tion. Therefore, an injector location does not affect the 
oscillation pattern as we observe. 

In our analysis of the snake states so far, we have re¬ 


stricted ourselves to a semiclassical picture. Since most 
of the oscillations we observe occur at high filling factors 
(up to 150, see Supplementary Fig. 5), the semiclassical 
picture is applicable. However, at low filling factors, evi¬ 
dence of the SdH oscillations, which requires a quantum 
approach, can be seen around (Vbg? Vtg) = (0, 0) V in 
Fig. 2c and Supplementary Fig. 3f. These SdH oscilla¬ 
tions disappear at high density and the oscillations due 
to the snake states become dominant. 


In a second device we measure, there is evidence that 
both edge states and snake states can coexist (see Supple¬ 
mentary Fig. 6). At low field, we observe the oscillations 
due to the snake states, similar to the first device. As we 
increase magnetic field, SdH oscillations and the oscilla¬ 
tions due to the snake states are both apparent and cross 
each other which imply the coexistence of edge states and 
snake states. Further study is required to understand the 
interplay between the edge states and the snake states 
along a p-n interface p32l 134] and also how energy gaps of 
the Landau levels affect the snake states [2111231 f35] . 
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Methods 


Figure lb shows an AFM image from one of our de¬ 
vices before transferring the top hBN. The position of 
the top gate is shown as an orange region. We first 
transfer graphene onto high-quality hBN [36] and subse¬ 
quently pattern it into a Hah bar shape in oxygen plasma. 
Chromium and gold are used for contacts and top gate, 
defined by electron-beam lithography. We anneal the 
sample in forming gas to reduce residues from fabrica¬ 
tion. In addition, for the device shown in Fig. lb, we 
also perform mechanical cleaning[37: to further clean the 
device. 

We measure conductance by applying a small voltage- 
bias excitation (< 2 mV) at 13 Hz to an injector and 
measuring the current through a collector. The mag¬ 
netic held is applied normal to graphene. Ah measure¬ 
ments were done at 5 K unless otherwise stated. A held 
effect mobility of our sample is r\j 40,000 cm 2 V 1 s h 
The capacitive couplings of back gate and top gate are 
Cbg = 9 x 10 -5 Fm“ 2 and Ctg = 3 x 10 -3 Fm“ 2 , 
determined from Shubnikov-de Haas oscillations using 4- 
probe measurement at high magnetic held. We observe 
the conductance oscillations due to snake states in two 
devices (see Supplementary Fig. 6 for data from a second 
device). 
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